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Abstract 

In this paper, we study minimization problems among Bravais lattices for finite energy per 
point. We prove - as claimed by Cohn and Kumar - that if a function is completely monotonic, 
then the triangular lattice minimizes energy per particle among Bravais lattices with density 
fixed for any density. Furthermore we give an example of convex decreasing positive potential 
for which triangular lattice is not a minimizer for some densities. We use the Montgomery 
method presented in our previous work to prove minimality of triangular lattice among Bravais 
lattices at high density for some general potentials. Finally, we deduce global minimality among 
all Bravais lattices, i.e. without a density constraint, of a triangular lattice for some parameters 
of Lennard-Jones type potentials and attractive-repulsive Yukawa potentials. 
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1 Introduction and statement of the main results 


The two-dimensional crystallization phenomenon - that is to say the formation of periodic struc¬ 
tures in matter, most of the time at very low temperatures, - is well known and observed. For 
instance, similarly to 45], the following may be mentioned : Langmuir monolayers, Wigner crystal, 
rare gas atoms adsorbed on graphite, colloidal suspensions, dusty plasma and, from another point 
of view, vortices in superconductors. In all these cases, particle interactions are complex (quantum 
effects, kinetic energy, forces related to the environment) and this implies that the physical and 
mathematical understanding of this kind of problem is highly complicated. However, we would like 
to know the precise mechanisms that favour the emergence of these periodic structures in order to 
predict crystal shapes or to build new materials. 


Semiempirical model potential with experimentally determined parameters are widely used in vari¬ 
ous physical and chemical problems, and for instance in Monte Carlo simulation studies of clusters 
and condensed matter. A widespread model is the radial potential, also called “two-body poten¬ 
tial”, which corresponds to interaction only depending on distances between particles. This kind 
of potential, based on approximations, seems to be effective to show the behaviour of matter at 
very low temperature, when potential energy dominates the others. There are many examples, 
that can be found in but they are usually constructed, except for very simple models such as 
Hard-sphere, with inverse power laws and exponential functions, easily calculated with a computer 
if we consider a very large number of particles. For instance we can cite : 


the Lennard-Jones potential r 


0-2 ai 1 , 

-, where the attractive term corresponds to the 

^X2 J-Xl ’ ^ 

dispersion dipole-dipole (van der Waals : ~ r“®) interaction, initially proposed by Lennard- 
Jones in 3^ to study the thermodynamic properties of rare gases and now widely used to 
study various systems, the best know being for {xi,X 2 ) = (6,12); 

the Buckingham potential r i-t- -1 proposed by Buckingham in [3 and including 

rpO 

attractive terms due to the dispersion dipole-dipole (~ r“®) and dipole-quadrupole (~ r“®) 
interactions, and repulsive terms approximated by an exponential function; 

g-br 

the purely repulsive screened Coulomb potential r i-t- a— 


, also called “Yukawa potential”, 

proposed by Bohr in for short atom-atom distances and used for describing interactions 
in colloidal suspensions, dusty plasmas and Thomas-Fermi model for solids [1, [^; 


^In a system of interacting electrons, where the coulomb interaction energy between them sufficiently dominates 
the kinetic energy or thermal fluctuations 
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the Born-Mayer potential r ^ ae used by Born and Mayer in in their study of the 
properties of ionic crystals in order to describe the repulsion of closed shells of ions. 


Many mathematical worked were conducted with various assumptions on particles interaction ; 
hard sphere potentials (dill . 48|; oscillating potentials 57]; radial (parametrized or not) potentials 
f^ . [^ . [l^ . Im ]: molecular siniulations with radial potentials B, HH, 0; three-body (radial and 

radial potentials and crystallization amoi^ Bravais lattices 


angle parts) potentials [371 . 


(Number Theory results and applications) S fill. B- 27. 17. fisl . 16, 54, It 


vortices, in 

superconductors, among Bravais lattice configurations [52l. l65l. l3| . Writing these problems in terms 
of energy minimization is common to all these studies. Furthermore, in many cases, triangular 
lattices (also called “Abrikosov lattices” in Ginzburg-Landau theory [l|, or sometimes “hexagonal 
lattices”), which achieves the best-packing configuration in two dimensions, is a minimizer for the 
corresponding energy. 



Fig. 1 : Triangular lattice 


A clue 
lattice 


to understanding this optimality, which is claimed in 14, p. ISojl, is the fact that triangular 
minimizes, among Bravais lattices, at fixed density, energies 


L^Ef[L\-.= Y. fiWpf) 

peL\{0} 


where ||.|| denotes Euclidean norm in and / : —)• M is a completely monotonic function, 

i.e. Vk G N,Vr G (—1)^ > 0. Moreover, Cohn and Kumar conjectured, in IJ, Conjec¬ 

ture 9.4], that the triangular lattice seems to minimize energies Ef among complex lattices, i.e. 
union of Bravais lattices, with a fixed density. Hence, it is not surprising, as for the Lennard-Jones 
potential we studied in 0, that some non-convex sums of completely monotonic functions give 
triangular minimizer for their energies at high fixed density. We observed this behaviour in works 
of Torquato et al. 3^, M]. However, it is important to distinguish mathematical results and phys¬ 
ical consistency. Indeed, at very high density i.e. when particles are sufficiently close, kinetic and 


^We cite only papers about 2D problems. 

proof of this assertion will be given in Section 3.1. 


3 
























quantum effects cannot be ignored and our model fails. For instance, Wigner crystal appears if the 
density is sufficiently low and matter obviously cannot be too condensed. Nevertheless, this kind 
of result is interesting, whether in Number Theory or in Mathematical Physics and this study of 
energy among Bravais lattices is the first important step in the search for global ground state, i.e. 
minimizer among all configurations. For instance, we have recently found in a deep connexion 
between behaviour of vortices in the Ginzburg-Landau theory, and more precisely works of Sandier 
and Serfaty [^, 53, [s^, and optimal logarithmic energy on the unit sphere related to Smale 7th 


Problem. Thus the optimality of triangular lattice, among Bravais lattices, for a renormalized 
energy W, which is a kind of Coulomb energy between points in the whole plane, gives important 
information about optimal asymptotic expansion of spherical logarithmic energy thanks to works 
by Saff et ah [i^, [^. 


The aim of this paper is to prove this minimality of triangular lattice at high density, with the 
same strategy as in our previous work [^, that is to say the use of Montgomery result [4^ about 
optimality of triangular lattice at a fixed density for theta functions 

L^0i(a):= ^ 

peAfo} 

for some general admissibl^ potentials /, summable on lattices and such that their inverse Laplace 
transform fif exists on [0,+oo). Hence, as in the classical “Riemann’s trick” that we used in 0], 
we can write an integral representation of energy Ef which we deduce a sufficient condition for 
minimality of triangular lattice among Bravais lattices of fixed densit}!!. This is precisely the aim 
of our first main theorem, which we now state. 


Theorem 1.1. For any admissible potential f, for any H > 0 and any Bravais lattice L of area A, 
there exists a constant Ca, which not depends on L, such that 


TT 


r-\-oo 


Ef[L] = j 


Ol 


2AJ 


y V/1 -T) + ^/ (^ 


TT 




yAj 


Try 


A 


dy + Ca 


( 1 . 1 ) 


where pLf is the inverse Laplace transform of f. Moreover, if 

(^) ^ (t) - ° 

/ 2A r 

then the triangular lattice of area A, i.e. Ka = \ —^ Z(l, 0) 0 Z(l/2, \/3/2) 

y vS L 

minimizer of L ca Ef[L\, up to rotation, among Bravais lattices of fixed area A. 


( 1 . 2 ) 

, is the unique 


Sufficient condition can be applied for some general functions /. More precisely we will 

consider the following potential^, defined for r > 0, which we will explain the interest throughout 
the paper : 

• Sums of screened coulombian potentials : 


Fa,x{r) = 

i=l 

rigorous definition will be given in preliminaries. 

® Actually, as in [J , we will write all our results in terms of area, that is to say the inverse of the density. 

®It is important to distinguish potential / and the function r i—^ f(r^) that we sum on lattices to compute Ef. 
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n 

with 0 < xi < X 2 --- < Xn, CLi £ K* for all 1 < i < n and a* > 0; 

i=l 


Sums of inverse power laws : 




2=1 


with 1 < xi < X 2 < ... < Xn, CLi £ IK* for all 1 < i < n and Un > 0; 


• Potentials with exponential decay : 

m 

2=1 

with 3/2 < xi < X 2 < ... < Xn, CLi £ IK^ for all 1 < i < n, a„ > 0, bj,tj £ M* for all 1 < j < m. 

Thus, even though our method is not optimal, we will give explicit area bounds in Propositions 15.ll 
16.3116.101 and 17.21 with respect to potential parameters, above which minimizer is triangular and 
we give conditions on parameters, for potentials tpn,x ^-rid in order to get a triangular global 
minimizer, i.e. without area constraint, in particular when the potential has a well. This is the aim 
of our second theorem. 


Theorem 1.2. Let functions ‘Pa,x, ,ya,x,yjfx fa,x,b,t be defined as before. 

A. Minimality at high density. If f & {pa,x,ya,x, fa,x,b,t} then there exists Aq > 0 such that for 
any 0 < A < Aq, Aa is the unique minimizer, up to rotation, of L among Bravais lattices 

of fixed area A. 


B. Global optimality without an area constraint. We have the following two cases 
1. Let be the attractive-repulsive potential defined by 

p-X2r p-xir 

AR/ \ ^ ^ 

Pa,x y) = a 2 — -«1— 


where 0 < oi < 02 and 0 < xi < X 2 . //ai, 02 , xi, X 2 satisfy 


--^eV ^2J > 1 

02(1 + vr) 


ai (aiX2 + XI (02 - ai)7r) > 1 

02X2 (ai + (02 - ai) 7 r) “ 


(1.3) 


then the minimizer of L i-a E^aAI'] among all Bravais lattices is unique, up to rotation, and 
triangular. In particular it is true if a 2 = 2ai and xi < 0.695x2. 


2. Let Vnx be the Lennard-Jones type potential defined by 



02 Oi 


with 1 < xi < X 2 and ( 01 , 02 ) £ (0,+oo)^. ITe set h{t) = 7r~^r{t)t. If h{x 2 ) < h{xi) then 
the minimizer La^x of L i-a EyLj[L] among all Bravais lattices is unique, up to rotation, and 
triangular. Moreover its area is 


_ / O2X2CA1 (2X2) \ “^ 2—1 
' “ Uia^iCAi(2xi)y' 

In particular, it is true if (xi, X 2 ) £ {(1.5, 2); (1.5, 2.5); (1.5,3); (2, 2.5); (2; 3)}[l|. 


^See Section 6.4 for numerical values. 
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We proceed as follows, we start below with some preliminaries where we recall Montgomery result 
about optimality of Ka for theta functions 9^ and we give the definition of admissible potential. 
Then we prove in Section 3 the optimality of for every A when / is completely monotonic and 
we give an example of strictly convex, decreasing and positive potential V such that Ka is not a 
minimizer oi Ef for some A. Theorem 11.11 is proved in Section 4, with some general applications. 
Furthermore we discuss optimality and improvement of this method. Finally we prove our Theorem 
fL2] in next sections where we present the interest, in molecular simulation, of studied potentials 
and we prove additional results. Throughout the paper, we give numerical results and examples. 


2 Preliminaries 

2.1 Bravais lattices, zeta functions and theta functions 

We briefly recall our notations in [^. Throughout this paper, ||.|| will denote the Euclidean norm in 
M^. Let L = ZmSZi! be a Bravais lattice of M^, then by Engel’s theorem (see 20]), we can choose 
u and V such that ||tt|| < ||n|| and (S7T) G 


to rotations and translations and the fact t 


3’ 2J 


in order to obtain the unicity of the lattice, up 


lat the lattice is parametrized by its both first lengths 


||n|| and ||t;||. We note |L| = jjn A ?;|| = ||n||||n|| |sin(u7S)| the aresu of L which is in fact the area 

/ 2A r 1 

of its primitive cell. Let = > —= Z(l, 0) 0 Z(l/2, y/2>l2) be the triangular lattice of area A, 

y \/3 1 J 

then ||tt|| is called the length of this lattice. 

For real s > 2, the Epstein zeta function of a Bravais lattice L is defined by 

p&L* 

where L* := L\{0}. As proved in [l^ . Proposition 10.5.5 and Proposition 10.5.7], we can write 
Cl{s) in term of L-function or Hurwitz zeta-function. More precisely, for L = 7? and L = Ai the 
triangular lattice of area 1, we have, for any s > 1, 


Cz2 (2s) = 4L_4 (s)C(s) = [C(5,1/4) - C(5,3/4)1 , 


( 2 . 1 ) 


Cm(2s) = 6(^) C(s)L_3(.) = 6(^) 3-C(^)[C(., 1/3)-C(., 2/3)1, (2-2) 


+ 00 


where C is the classical Riemann zeta function ^(s) := Ld defined by 


2=1 


H-cxd 


Ld{s) ■■= ^ 


n=l 


D 


n 


n 


is the Dirichlet L-function associated to quadratic field Q{iy/—D), with ^— ) the Legendre symbol, 
and, for x > 0, 

+CXD 


D 


C(s,x) := ^(n + : 


n=0 


is the Hurwitz zeta function. Hence both these special values are easily computable. 


“We choose, as in to write results in terms of area and not in terms of density (which is its inverse). 
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Now we recall fundamental Montgomery’s Theorem about optimality of Ka among Bravais lat¬ 
tices for theta functions : 

Theorem 2.1. (Montgomery, [i^ ]! For any real number a > 0 and a Bravais lattice L, let 


OUa) :=0i(^a) = ^e-2™M^ 

p&L 


(2.3) 


where Ql is the Jacobi theta function of the lattice L defined for Im{z) > 0. Then, for any 
a > 0, h-A is the unique minimizer of L ^ 0L{a), up to rotation, among Bravais lattiees of area A. 


Remark 2.2. This result implies that the triangular lattice is the unique minimizer, up to rotation, 
of L I— )• Cl{s) among Bravais lattices with density fixed for any s > 2 which is also proved by Rankin 
in [^, Cassels in 111, Ennola i n and Diananda in 181. Montgomery deduced this fact by the 
famous “Riemann’s trick” (see jsjor for a proof): for any L such that D = 1, 


for Re(s) > 1, CL(2s)r(s)(27r) ^ = 


^ - 1)(«^ + «'-^) —. (2.4) 

s — 1 s ./i a 


In Section 4.2 we will prove general Riemann’s trick 11.11 which we call integral representation of 
energy, for admissible potentials in order to use Montgomery method in a general case. 


Remark 2.3. We can find in |44l . Appendice A] other proof of minimality of some theta functions 
based on result of Osgood, Phillips and Sarnak [^, Corollary 1(b) and Section 4] about Laplacian’s 
determinant of flat torus, which has some deep connection with other energies (for instance, see [^, 
Theorem 2.3]). 


2.2 Admissible potential, inverse Laplace transform and lattice energies 
Definition 2.1. We say that f : {Re{z) > 0} —?■ M zs admissible if : 

1. there exists rj > 1 sueh that \f{z)\ = 0{\z\~^) as \z\ -too; 

2. f is analytie on {z € C] Re{z) > 0}. 

If f is admissible, we define, for any Bravais lattiee L ofE?, 


Efir := /(llpit) 

pGL* 

which is the quadratic energy per point of lattice L created by potential f . 


Remark 2.4. As a consequence of 47|, Theorem 5.17, Theorem 5.18], we get, by direct application 
of inversion integral formula : 

• There exists an unique inverse Laplace transform //Jl, which is continuous on (0, -too); 

• We have /Uj(0) = 0. 

Remark 2.5. This definition excludes two-dimensional Coulomb potential r i—)• — logr because all 
its quadratic energies are infinite. However we can define a renormalized energy as in 5^ or in 


^We will sometimes write C and C ^ for Laplace and inverse Laplace operators. 
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2.3 Completely monotonic functions 

The class of completely monotonic functions is central in our work. Indeed, as we will see in Sec. 
3, these functions have good properties for our problem of minimization among lattices with fixed 
area thanks to the Montgomery theorem 12.11 


Definition 2.2. A function f : (0,+oo) —)• M_|_ is said to be completely monotonic if, for 
any k gN and any r > 0, 




Examples 2.6. We can find a lot of examples of completely monotonic functions in Here we 
give only some interesting classical admissible potentials / : 


• 14 (r) = r X > 1; 


n 

• U,.(r) = a^r where a* > 0 and Xj > 1 for all i; 
i=l 


fa{r) = e , a > 0, a G (0,1], see Corollary 1]; 


• Modified Bessel function, i.e. one of the two solutions of r'^y” + ry' — (r^ + i'‘^)y = 0 which 
goes to 0 at infinity, is Ky{r) = cosh(nt)dt, n G R. Moreover, r i-)- Ky{y/r) is 

also completely monotonic (despite we thought in [J]). 


g-av/r 

• Vsc{r) = —a > 0; 
vT 


• ^a{r) 



a > 0. 


Remark 2.7. We remark that if r i—)• /(r) is completely monotonic, it is not generally the case 
for r I—)• /(r^). For instance r i—>■ is completely monotonic, but r does not check this 

property. 


Now we give the famous connection between completely monotonic function and Laplace transform 
due to Bernstein in [^. 

Theorem 2.8. (Hausdorff-Bernstein-Widder Theorem) A function f : R is com¬ 

pletely monotonic on M+ if and only if it is the Laplace transform of a finite non-negative Borel 
measure /i on M+, i.e. 

r+oo 

fir) = CMr)= / e-^^dpi{t). 

Jo 

Remark 2.9. If / is admissible and completely monotonic, then 


dy{t) = p,f{t)dt and > 0, a.e. on (0,+oo). 

Remark 2.10. Actually Schoenberg proved in that r i—)• /(r) is completely monotonic if and 

only if r I— >• /(r^) is a positive definite function in M, i.e. for any N G N\{0,1}, any xi, ...,xn G M 
and any ci,..., cm G M, we have 

N 

CiCjf{\xi - Xjp) > 0 

or, by Bochner Theorem (see j^), if and only if r i—>■ /(r^) is the Fourier transform of a positive 
finite Borel measure on R. 






Positivity of Fourier transform of a radial potential is a key point in crystallization problems. 
Indeed Nijboer and Ventevogel proved in that it is a necessary condition for a periodic ground 
state (Bravais lattices) and Siito studied in his work potentials / such that f{k) > 0 and 
f{k) = 0 for any ||fe|| > Rq and proved some interesting crystallization results at high densities. 
Unfortunately, as Likos explained in [^, this kind of potential, oscillating and with inverse power 
law decay, seems to be difficult to achieve physically. 


Actually it is more common to use Fourier transform in problems of minimization of lattice energy 
because we have the Poisson summation formula and the natural periodicity of sinus and cosinus. 
Furthermore, applications of classical formula allows to obtain some interesting results, as in IJ, 


Proposition 9.3]. However we will show in Section 4 that inverse Laplace transform also seems well 
adapted to our problem and gives simple calculations. Indeed, Fourier methods as in 14, 5^ is 

good for more general minimization problems and our method is a better choice for minimization 
among Bravais lattices because of integral representation m- 


2.4 Cauchy’s bound for positive root of a polynomial 

In this part, we recall Cauchy’s rule explained in Note III, Scolie 3, page 388] for upper bound 
of polynomial’s positive roots (see also for simple proof). 

Theorem 2.11. (Cauchy’s rule) Let P a polynomial of degree n> 0 defined by 


P{X) = Y,ociX\ a„>0 

j=0 


where a, < 0 for at least one i, 0 < i < n — 1. If X is the number of negative coefficients, then an 
upper bound on the values of the positive roots of P is given by 


Mp = max 


-Xotj 

CXn 


1 

n — i 


Remark 2.12. This Theorem stays true for upper bound on the values of the positive zero of any 
function p defined by 

n 

p{y) = '^Oiiy''\ > 0 

i=l 

where 0 < < ... < t'n are real numbers and we obtain 


Mp = max 

v,ai<0 



(2.5) 


This result will be useful for technical reasons in the following sections, because we will want 
positive zeros less than 1 to apply our sufficient condition inTheorem 11.11 and to prove Theorem 

OA. 


3 Completely monotonic functions and optimality of 

In this part we begin to state a simple fact connecting positivity of inverse Laplace transform and 
minimality among lattices at fixed area. Furthermore we will give an example of strictly convex, 
decreasing, positive potential for which there exists areas so that the triangular lattice is not a 
minimizer among Bravais lattices with fixed area. 
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3.1 Optimality at any density 

The following proposition, claimed by Cohn and Kumar in Q, page 139], is a natural consequence 
of Montgomery and Hausdorff-Bernstein-Widder Theorems. 

Proposition 3.1. Let f be an admissible potential. If f is completely monotonic then, for any 
A> Q, Ka is the unique minimizer, up to rotation, of 


LeAEf[L]= Y^fm' 

p&L* 


among lattiees of fixed area A. 

Proof. As / is admissible, we can write, 


f + OO 

f{r) = / 

Jo 


r. — Tt 


Pf{t)dt 


and it follows that 


E/W = E /dipii') = E / 

r+oo 

Jo 


o- 4 p\\ 


, ^+oo 

'hf{t)dt = Y 
Jo peL* 


-t\\p\\‘ 


pLf{f)dt 


Pf{t)dt 


and 


r+oo 

Ef[L]-Ef[AA] = 

Jo 


9l 


27r 


- 


t 


Hf{t)dt. 


If / is completely monotonic, by Theorem 12.81 yLf{r) > 0 for almost every r E (0, +00). Moreover, 
by Montgomery Theorem 12.11 for any t > 0 and any Bravais lattice L of area A, 

9l ( 7 ;-^ - > 0, 


27r 


27r 


hence Ef[L] > Ef[AA] for any L such that |L| = A and is the unique minimizer of the energy 
among Bravais lattices of fixed area A. □ 


Remark 3.2. We can imagine that the reciprocal is true, i.e. if / is not completely monotonic, then 
there exists Aq such that Ayi^ is not a minimizer among Bravais lattices of area Aq fixed. In next 
subsection will give an explicit example correlated with Marcotte, Stillinger and Torquato results 
in about the existence of unusual ground states with convex decreasing positive potential. 

Examples 3.3. A direct consequence of this theorem is the minimality of triangular lattice among 
lattices for any fixed area for the following energies : 

• Ev„,[L] = X > 1 is the first natural example given by Montgomery in (43l |. 

n 

• Ev„[ti = E aiCL{‘Jxi) where Oj > 0 and Xj > 1 for all i, 

i=l 

• Ef^[L] = Y^ E (0,1], in particular Ef^^^[L] = Y^ 

pSL* pSL* 

• 2 ATiydlpll) , zz E M which generalizes our study of lattice energy with potential 

peL* 

Kq in in Thomas-Fermi model case; 


10 






o > 0, which corresponds to lattice 


potential interaction and can explain formation of triangular 


26 i; 


energy for screened Coulomb 
Wigner crystal at low density 


• E^AL] = Y. 

p&L* 



a > 0. 


3.2 Repulsive potential and triangular lattice 

In this section we give an example of stricly convex decreasing positive radial potential V so that, 
for some densities, a minimizer of Ey among Bravais lattices of density fixed cannot be triangular. 
As Ventevogel and Nijboer proved in 6l|, a convex decreasing positive potential allows to obtain, 
in one dimension and for any fixed density, a dilated of lattice Z as unique minimizer among all 
configurations. Thus the two-dimensional case is deeply different. 


Let 


be the potential and 


V{t) 


14 



Ey[L\ = 14a (4) - 40a(6) + 35a (8) 


the energy per point of a Bravais lattice L. 


(3.1) 



Proposition 3.4. (Strictly convex potential and non optimality of triangular lattice) Let 

V be given by then 

• V is strictly positive, strictly decreasing and strictly convex on (0,-|-oo); 

• There exists Ai, A 2 such that A a is not a minimizer of Ey among all Bravais lattices of 
area A G (^ 1 ,^ 2 ). 

Proof. We have 

,,, , 14r2-40r + 35 
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and the discriminant of polynomial — 40X + 35 is Ai = —360 < 0, hence V{r) >0 on (0, +oo). 

We compute 

^ -4(7r^-30r + 35) 

and the discriminant of — 30A + 35 is A 2 = —80 < 0, therefore V'{r) < 0, i.e. V is strictly 
decreasing on ( 0 ,+oo). 

Moreover, we have 

, 4(21r2 - 120r + 175) 

W =- ;:6 - 

and the discriminant of 21A^ — 120A + 175 is A 3 = —300 < 0, then V"(r) > 0 on (0, + 00 ), i.e. V 
is strictly convex on ( 0 ,+ 00 ). 

For the second point, we have the following equivalences 


Ev[L] > for any \L\ = A 

^ 14a(4) - 40 Cl(6) + 35 a(8) > 14 Ca^(4) - 40Ca^(6) + 35Ca^( 8) > 0 for any \L\ = A 
^ “ Cai(4)) + ^ (Cl(6) - Cai(6)) + ^ (Cl(8) - Cai( 8)) > 0 for any \L\ = 1 

^ 14(Cl(4)-Cai(4))A2-40(Cl(6)-Cai(6))71 + 35(Cl(8)-Cai( 8)) >0for any |L| = 1 
Pl{A) > 0 for any \L\ = 1 


where the discriminant of Pl(A) = 14(Cl(4) - Cai (4)) -40 (Cl (6) - Cai (6)) A+35 (Cl(8) -Cai(8)) 

is 

A(L) = 1600 (Cl(6) - Cai(6))2 - I960 (Cl (4) - Cm (4)) (Cl ( 8 ) - Cm ( 8 )). 


For L = 1? the square lattice of area 1, we obtain A(Z^) ~ 24.231435 > 0 then there exist two 
positive numbers Ai and yl 2 such that Pj 2 {A) < 0 for any Ai < A < A 2 . Hence, is not a 
minimizer of Ey among Bravais lattices with fixed area A if Ai < A < A 2 . More precisely we get 


Ai PS 2.3152307 and A 2 Ri 3.759353. 


□ 

Remark 3.5. It follows, from the previous proof, that function r 1 —)■ V{r‘^) is also strictly positive, 
strictly decreasing and strictly convex on ( 0 ,+oo). 

Remark 3.6. Actually, the previous proof implies that, for any A G (Ai, A 2 ), 

EvIVAZ^] < Ev[Aa]. 

Moreover, this interval seems numerically to be optimal, i.e. for any A 0 [Ai, A 2 ], A a seems to be 
the unique minimizer, up to rotation, of L 1 —)■ Ef[L] among Bravais lattices of fixed area A. 

4 Sufficient condition and first applications 

Now we study the case of non completely monotonic potential /, i.e. ///is negative on a subset of 
(0, + 00 ) of positive Lebesgue measure. 
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4.1 Integral representation and sufficient condition : Proof of Theorem 11.11 

Proof. Let L be a Bravais lattice of area A and / an admissible potential. Firstly we prove the 
integral representation [TTT] of energy Ef[L] : 


Ef[L] := ^ /(Iblp 

p&L* 


r 


9l 

r*+00 


JL 

2A 




Try 




dy 


Indeed, for a Bravais lattice L of with \L\ = 1/2, we have, as in 0, by t = 2?™, u = y ^ and 
Montgomery’s identity = ydLiu) (proved in [i^) : 


/ + CO ^ + CXD 

p^n- p^n- p&L* 

r+oo 

/ [6l{u) — 1] yLf{2Tru)du 

Jo 

/ I r+oo 

[0L(tt) — 1]/rj(27rn)(in + 27r J [6l{u) — 1] fif{27ru)du 

~ ^ [9l{u)- l]y,f{2-jru)du 

Mi(y) “ 1 ]/^/^ + - l]^/(27rM)du 

^ [0L(?r) - l]^/(27rn)dn-27r^ M/(y) ^ 

^ [6'l( 2/) - 1] (^y"V/+/^/(2vry)^ dy+ 27r^ yf(^^'^{y~^-y~‘^)dy 


= 27r 

= 27r 

= 27r 

= 27T 

= 27r 

= 27T 


y)^ 

Now we have, by change of variable t = y~^, 

r* + OD 


i: 


—'] {y ^ -y ‘^)dy 
y J 


< / \yLf{2'Kt)\ {t — l)dt < +00 


because /.i/ is continuous on M^, /rj(0) = 0 and 1t ^ is integrable at the neighbourhood of 0. 
Hence, for L such that \L\ = A, we have 

EflL] = E /(Ill’ll") = E /(2-4||Pll") 


pGL* 


p&L* 


where L = V2AL, \L\ = 1 / 2 . By identities y/(fc.) = and 9j^{s) = 9l 


TT /■+°° 


9l 


JL 

2 d 


- 1 


y V/ 


—) 

yAj 


+ l^f{^))dy 


+ 


7T 


r+oo / 

f 


TT 

yA 


(y - y )dy 


and Ca ^ j yf ^ ~ y is a finite constant which not depends on L. Now our 

sufficient condition is clear because, for any Bravais lattice L of area A, we have 


/*+oo p 

Ef[L]-Ef[AA] = j 


dLirrj] -9 aa 


2A 


2AJi 


9A{y)dy. 
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By Montgomery theorem, Ol{u) — 9\^{u) > 0 for any n > 1 and any L. Thus, if 

(^) ^ (t) - ^ ^ ^ 

then it follows that 

Ef [L] — Ef [Ayi] > 0 

and is the unique minimizer of Ef, up to rotation, among lattices of fixed area A. □ 

4.2 Minimization at high density for differentiable inverse Laplace transform 

In this part we give two results, in the case of differentiable inverse Laplace transform, which are 
based on our Theorem 0 

Proposition 4.1. Let f be an admissible potential such that pLf is with derivative fr'f. If 
L Lf{y) > 0 on ^,+oo^, 

then A^ is the unique minimizer, of Ef, up to rotation, among Bravais lattices of fixed area A. 
Proof. We write, for any y > 1, 


with 


Therefore, we get 


<^a(v) » (^) + W / (^) . 


(T^y\ 

, Try 

' ! \ -3 t { 

V A J 

+ A 



Assumption 1. implies that pj > 0 for any y > 1. Moreover, it is clear that point 2. means 
that p'f {^y) — y~^y'f (i^) — 0 y ^ 1; hence u'j^{y) > 0 for any y >1. As 

ua{ 1) = 2p,f Q) > 0 

we have u^iy) > 0 for any y > 1 and it follows that 

gA{y) >0 Vy > 1 

and by Theorem 11.11 A^ is the unique minimizer, up to rotation, of Ef among Bravais lattices of 


hxed area A. 


□ 


Corollary 4.2. Iff is an admissible potential such that its inverse Laplace transform fj,f is convex 
on (0, +oo), then there exists Aq > 0 such that for any A E (0, Aq), Aa is the unique minimizer of 
Ef, up to rotation, among Bravais lattices of fixed area A. 

Proof. As /xj- is convex, there exists tq > 0 such that, for any r > tq, /U/(r) > 0. Moreover, for any 

1‘i (I'i') ^ '■/ {£} 

TT '^y 

because < — and /x/ is convex. Hence, as y~^ < 1 for any y > 1, we get both points 1. and 

7f 

2. of Proposition 14.II for any A such that 0 < A < Aq := —. □ 

ro 
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4.3 Remarks about our method 

As we saw in j^, for Lennard-Jones case, our method is not optimal to finding all areas such that 
is the unique minimizer, up to rotation, of Ef among Bravais lattices of fixed area A. The 
general and difficult main problem is to find all A such that, for any Bravais lattice L of area A, 


Ej[L]-Ef[AA] = jl^ 


vr /+-r„ (JL 


Ol 


\2A 


\2AJ\ 


gA{y)dy > 0 


where gA{y) '■= y d-f (”^)' imagine that even if qa is not positive almost 

everywhere on [l,+oo), the positive part of this integral can compensate the negative one. For 
instance, if we consider, as in [^, /(r) = — 2r~^, then 


gA{y) = 


TT 

A 2 


TT" 


A35! 


(/ + ?/■")-r-y 


-2 


TT 


and we plot graphs of y e-)- „ (y^ + y — y^ — y ^ for A = 0.8 (on the left) and A = 1 (on the 

A'^S! 

right). 


0 0.2 0.4 0.6 



1.4 1.6 1.8 2 



/ VA 

gA{y)dy where yA is the second zero of qa-, A € {0.8,1} 
Thus a fine study, with respect to lattices L and real y, of the behaviour of positive function 


y 


y 


is necessary. However we find it difficult at this time. Indeed, 
• for any Bravais lattice L of area A, 


lim /S.L{y) = 0 

y^+oo 


and Ai exponentially decreases ; 

• if complete monotonicity is a necessary condition to optimality of A^ for any fixed area A, 
then function y i-t- ^.^{y) is not decreasing on [1, +oo) for any A and any L with area A. 
Indeed, is decreasing on (1, +oo) if and only if for any t > 1, A'^{y) < 0, i.e. 

\/A,\/L,Al decreases on (l,+oo) 
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pSL* 


A 


|2 —2L 

e -4 


< 0 


peA* 


VA,VL,Vy>l, ^ 

peA^ 


|2g-ip||p||^ < ^ 
pGL* 


|2 —2L 

e ^ 


which would be not possible, because r ^ re is never completely monotonic for y > 1. 

Hence comparing lA.L{y)gA{y)dy and AL{y)gA{y)dy seems difficult, even improving our 
method is possible (see for numerical values). 


5 Sums of screened Coulomb potentials 

In this part, we give the first simple example of application of Theorem 11.11 We consider non 
convex sums of screened Coulomb potentials and we prove minimality of at high density and 
global minimality of a triangular lattice among all Bravais lattices. 

5.1 Definition and proof of Theorem I1.2L A for ^pa,x 

n 

Definition 5.1. Let n G N*. For coefficients a = (oi, ..., On) € such that a* > 0 and for 

i=\ 

x = {xi,...,Xn) G (M+)"', we define 


n 

Fa,x{r) ■■= 

i=l 


g-Xir 


r 


and we set tCn '■= 


< 0 >. 

. i=l J 


As a proof of Theorem I1.2I A for this potential, we purpose to give an explicit bound for the 
minimality at high density as follows 

Proposition 5.1. Assume 0 < xi < ... < and let A such that 

E '' '' ' n \ '^ ^ 

i=l I 


A < min < min 

k&Ka 


TT 


TT 

Xn. 


then A^ is the unique minimizer of up to rotation, among Bravais lattices of fixed area A. 

Proof. We compute easily, because C~^[r~^e~^^'^]{y) = l[a;j,+oo)(y) any x* > 0 and any y > 0, 


Z=1 


and it follows that, for any y > 1, 

1 

9A{y) ■ ~ y Ifxj. + Cio) 


yT=i 


(s 




+oo) 


i=l 


A 


^ n n 


16 







TT 

As, by assumption, A < —, we have, for any 1 < i < n — 1, 

TT TT 

> -- > 1 . 


Axi Axi+i 


Hence we get 


QAiy) = < 


(l + y ^)^ai if 1 <y < ^ 


i=l 

n 




i=l ^ i=l 
di 


E' 

k i=l 


if y > 


As Er=i «* > 0 and, for any k ^ JCa, X]i=i — 0 ) obtain 


Vy G 


1 , 


TT 

AXr, 


u( 


TT TT 


k^iCa 

n 


u 


TT 


,+oo K yA(y)> 0 . 


Now if A: G )Ca, as H < min ■. . , 

fce/Ca I^Xfe+i y Ya=i<^' 


\Axk+i Axk\ vAxi 

we get, for any y G ( ^ 


E7 + E“-£^E“- + E“<£« 


y 

i=l " j=l 


i=l i=l 


and it follows that gA{y) > 0 for any y > 1. By Theorem ll.il is the unique minimizer of 
up to rotation, among Bravais lattices of fixed area A. □ 


5.2 Global minimality : Proof of Theorem II. 2L B.1 

Now we focus on particular “attractive-repulsive” case 

• a = (— 01 , 02 ) where 0 < oi < 02 ; 

• X = (X 1 ,X 2 ) with 0 < Xl < X 2 . 

Therefore we define, for any y > 0, 

p-X2r p-x\r 

Af?/\ G C 

^a,xir) 1 = 02 ^— -Oi^—. 

Now, let us prove Theorem II. 2I B.1. 

Proof. Firstly we study variations of (px,a to prove the existence of global minimizer La^x among all 
Bravais lattices and upper bound Oa,® for its area. Afterward we prove that inequalities m are 
equivalent with 

aa,x < min j—, — (— - l) | • 

[X2 X2 \ai J j 

Thus, by direct application of Theorem 15.11 if A < min \ — 1) >, A^ is the unique mini- 

I 312 ^2 V ^1 J J 

mizer among Bravais lattices of fixed area A, therefore La,x is triangular and unique. 
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STEP 1 : Variations of function 

We have, for any r > 0, 

<x(r) = ^ [ai(l + - 02(1 + X 2 r)e-"^^] 

and it follows that 

> 0 ga,x{r) ■■= {X2 - xi)r + ln(l + Xir) - ln(l + X 2 r) + In > 0. 

As, for any r > 0, 

, ,, {X2 - xi) (xiX2r'^ + (xi + X2)r) 

=- (l+r..r)(l + ..0 -> "■ 

ga,x is an increasing function on (0, + 00 ). We have 02 > ai, therefore In < 0 and there exists 
aa,x such that 

Vr G (0,aa,a;],5a,a;(’’) < 0, and Vr > aa,x,fl'a,x(?’) > 0. 

Thus we get ipa,x is a decreasing function on (0, aa,x] and an increasing function on {ua^x-, + 00 ). 


STEP 2 : The existence of global minimizer for ^ 

Variations of function ipa x and the fact that lim ipaxix) = +00 and goes to 0 at infinity implies that 

’ r —>-0 ’ 

r >0 

global minimizer exists. Indeed, this problem can be viewed like a minimization problem of a three 
variables function. By previous limits we can restrict this problem with variables in a compact set, 
and by continuity this problem has a solution La^x- 


STEP 3 : Upper bound for \La^x\ and conclusion 

Let = '^Ua^x © '^Va^x- If ||^a,x|| > then a contraction of all distances yields a new lattice 

with smaller energy because, by STEP 1, r i-G ipa,x{x‘^) is an increasing function on (y/o^,+ 00 ). 
Moreover, if ||ua,x|| > then a contraction of also gives a lattice with less energy. Thus 

we have ||rta,a;|| < ||w,a;|| < ^OLa.x- Now, because \La,x\ < ||wa,a:||||w,a:||, we ge10 


I La,a: I ^ 0:a,x- 


Now it is not difficult to check that 


’'Pa ,01 


—) > 0 ^ 

X 2 J 02 ( 1 +7rj 


- > 1 


and 


’ V ^2 V «1 // 02X2 (ai + (02 - ai)7r) 


hence ()5.2p holds and La,x is unique and triangular by Theorem 15.11 as explained at the beginning 
of the proof. 


^°This argument appears in [3, Proposition 4.1, ii)] and in in order to prove that the distance between two 
animals in a swarm is less than a specific “confort distance” between them, which minimizes a certain function. 
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STEP 4 : Example 

If we take 02 = 2ai then 


ai (aiX2 + xi(a2 - ai)7r) _ 

02X2 (ai + (02 - ai) 7 r) 02(1+ 7 r) 2 ( 1 + 7 r)'''^ ' X2 




Xi , 1-^W 

(1 -\ -7r)e'^ "“2/ , 


ll±^ -X+TT > ^ 

2 ( 1 + 7r)^ 

fi((X) := -X + log(l + X) - log(2 + 27r) + TT > 0. 


Now we set X = —vr and our condition becomes -(-e > 1, which is equivalent with 

2(1 + vr) 


As g'{X) = > 0 on M+, then g decreases and there exists X > 0 such that g{X) = 0. 


□ 


Numerically, we found X > 2.186, hence if A < 2.186, which corresponds to —tt < 2.186, i.e 

X2 

2 

xi < - X 2 ~ 0.695825x2, then g{X) > 0. In particular, it is true if xi < 0.695x2. 

vr 

Example 5.2. For instance, we can choose (xi,X 2 ) = (1,2). Thus, global minimizer of 

o-m? 


L ^ ^ ^aAWpf) = 2 ai E 


ai E 


is unique, up to rotation, and triangular. Hence we can construct potential with arbitrary deep 
well (using parameter ai) and with triangular global minimizer. 



Remark 5.3. This kind of potential seems not to be used in molecular simulation but this predic¬ 
tion of triangular ground state could be observed in the future. Furthermore our Theorem ll.2[ B.l 
allows to better understand ground state for parametrized potential with repulsion at short distance 
and quick decay at large distance, as in [6fll | where Theil proved global minimality of a triangular 
lattice among all configurations if the potential’s well is sufficiently narrow, i.e. with repulsion and 
decay sufficiently strong. 
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6 Nonconvex sums of inverse power laws 

In this part, we generalize our result in [^, which tackled only the classical Lennard-Jones case, for 
any nonconvex sums of inverse power potentials, that is to say optimality of triangular lattice 
at for high densities and non-optimality of this one for low densities. Furthermore we show that 
our method allows to obtain global minimizer, i.e. minimizer among all Bravais lattices without 
constraint of area, of Lennard-Jones type energies with small parameters. 


6.1 Definition and proof of Theorem I1.2L A for 

Definition 6.1. Let n > 1 be an integer and, for a = (ai, ■..,an) £ (M*)” such that an > 0, and 
X = {xi,...,Xn) £ (IK-i-)”' such that 1 < xi < ... < Xn, let 


VaAr) = Y. 

i=l 


di 

rpX{ 


We set I- := {i;ai < 0}, /+ := {z;ai > 0} and Oi 
(otherwise 14 ,x is completely monotonic). 


——. Moreover we assume that I_ / 0 
r(xj) 


Remark 6.1. In order to minimize among lattices, we should assume > 0 because we 

have I4,x(r) ~ anr~^^ as r — 0. Indeed we have I4,x(r) —>■ J-oo as r — 0 and I4,a;(r) —>■ 0 as 
r ^ J-oo therefore there exists minimizer of Ey^ ^ among Bravais lattices with fixed area. If o„ < 0, 
it is sufficient to do ||u|| —)• 0 to get Ey^^^[L] ——oo. 


Example 6.2. This kind of potential is widely used in molecular simulation. Indeed, besides 
Lennard-Jones potentials that we will study in the next subsection, it is sometimes necessary to 
consider some modifications of it. For instance, the (12 — 6 — 4) potential proposed by Mason and 
Schamp in [^, defined by 

V(r) = 

describes the interaction of ions with neutral systems. For instance, in fullerene Cgo, this potential 
describes interaction between a carbon atom in the polyatomic ion and a buffer gas helium atom. 
An other example, proposed by Klein and Hanley in [3i 


28( 1 for description of rare gases, more 


precise than Lennard-Jones, is the potential defined, for m > 8, by 


V(r) = ^ ^ ^ 

As in previous section, we give an explicit bounds for the minimality of A^ at high density in the 
following proposition. 


Proposition 6.3. If we have 


A < Trmin 
i&I- 




2tj{/_}|ai|F(a:^ 


xn — x^ 


( 6 . 1 ) 


then A A is the unique minimizer of Ey^^^, up to rotation, among Bravais lattices of fixed area A. 
Proof. By usual formula, we have 


hVaAv) 


E 

i=l 


yXi-l 


r(x. 
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and it follows that 


gA{y) := y Vi4,. + yv.,. (^) = jz ') 


2=1 




We set 

n 

PaAv) ■■= E 

2=1 

We notice that the term of high order is An-i with a„ > 0 and the number of negative 

coefficients is 2jj{/_}. Thus, by Cauchy’s rule 12.111 and more precisely generalization [23I an upper 
bound on the values of the positive zero of pa^x is 


Mp,,. ■= max ■ 

l£l- 


(2^{lA\ai\A^ 


O^r. 






dr; 


Xn + XA — 1 


because 2Xn — 1 — (Xn — Xi) = Xn + Xi — 1 and 2Xn — 1 — (Xn + Xi — 1) = Xn — Xi- 
We notice that 


A < 


TT mm 


(^n^(Xi 


iei- I V2tJ{/-}|ai|r(xn) 

dn 


xn — x^ 


= mm 
iei- 


dri 




djT, 

Mp^,^ < 1 


, Vi G /_ 
< 1, Vi G I- 


therefore the assumption implies that the largest zero of pa,x is less than 1. As > 0, it follows 
that Pa,x{y) — ^ any y > and then gA{y) > 0 for any y >l and by Theorem 11.11 if (|6.ip 

holds, then Au is the unique minimizer of Ey^ ^ among Bravais lattices of fixed area A. □ 


Remark 6.4. This result seems to be natural because for r close to 0, Va^xix) ~ anr~^'^ and for any 
A, Ka is the unique minimizer of L i—)• (^L{2xn) among Bravais lattices of fixed area A. However, if 
we fix A, ||m|| and ||u|| can be as larger as we want and the behavior of Va^x can be unusual. 
Furthermore, in the case 




Ql 

rpX\ 


where 01,03 are positive and 02 negative, our bound (16. ip does not depend on oi. For instance, if 
a = {p, —3,1) and x = (2,4,6), then, for any p, 


TT min 
i&I- 


anT{x^i) /r(4)y/2 

tt{i_}|oi|r(x.); / '^Ur(6)y' 


0.2867869 


which corresponds to triangular lattices of length ~ 0.5754589. 
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Example 6.5. For our counterexample (I.S.lh . i.e. V{r) = -^ H—j, a = (14,-40,35), 

X = (2,3,4) and ^{I-} = 1, hence 


vr min 
i&I- 


( anr(xi) /35r(3)y 

U{/-}|a.|r(x.); /“’^Uor(4); 


Ttt 

48 


0.4581488, 


which corresponds to triangular lattice of length ps 0.7273408. Thus, for A < —, is the unique 

48 

minimizer of Ey, up to rotation, among Bravais lattices of fixed area A. 

6.2 Non-optimality of at low density 
Proposition 6.6. If ai <0 and 


A > inf max 

^^■^1 2E/4- 
\L\ = 1 


f ll{-^+}«i(CL(2xi) - CAi(2xi)) \ 
V I«i|(Cl(22;i) - Cai(2xi)) ) 


( 6 . 2 ) 


that is to say if A is sufficiently large, then Aa is not a minimizer of among Bravais lattices 
of fixed area A. 

Proof. Let La = ^/ALi be a Bravais lattice of area A, with |Li| = 1, then 


n 

kVa.AJ^A] - -E^Va,r.[^A] = ^ Oj(C aa (2^:* ) “ Ci,A(2a^*)) 

i=\ 

n 

n 

= A-""^a,(CAi(2x,) -ai( 2x,))A—«. 
i=l 


We set 

n 

Pa,x,Lii^) ■= '^ai{CAii‘IXi) - CLi(2Xi))A^""^^ 

i=l 

As oi < 0 and, for any s > 1, 

Cai( 2 s)-Cli( 2 s) < 0 , 
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because Ai is the unique minimizer of L i—C l( 2s) among Bravais lattices of area 1, we can apply 
Cauchy’s rule 12.111 and more precisely its generalization (12.51) . The number of negative coefficient 
of Pa,x,Li is exactly and an upper bound on the values of the positive zero of Pa,x,Li for given 


Li, is 


Mp,,ALi) ■■= max 
2 E/+ 


f tl{-^+}Q»(CLi(2xt) - CAi(2xi)) \ i 

V |ai|(CLi(2a;i) - Cai(2xi)) J 


Hence, for any L such that \L\ = 1, if H > Mp^^^{L) then Pa,x,L{^) > 0. We conclude that if 
(j6.2p holds, then Ev^^^[^a\ — > 0 and cannot be a minimizer of among Bravais 

lattices of fixed area A. □ 


Remark 6.7. To compute explicitly a lower bound for A such that A^ is not a minimizer of 
energy Ey^^^, we can take L = in (16.21) and use equalities (12.11) and (|2.2I) (see next subsection 
for computations in Lennard-Jones case). 


6.3 Lennard-Jones type potentials : proofs of Theorems I1.2L A and I1.2L B.2 and 
numerical results 

Now we want to study more precisely the class of Lennard-Jones type potential. In we studied 
classical (12 — 6) Lennard-Jones potential Vlj(x) = — 2r“®, such that its minimizer is 1, and 

we proved that the minimizer of its energy among lattices with fixed area A is triangular for small 
A and it cannot be triangular for large A. Here we prove that our method gives interesting results 
for this kind of potential. 


Let 1 < xi < X 2 and ai,a 2 G (0, -|-oo), we define Lennard-Jones type potential by 

V^iir) := 


- “2 “1 


'pX‘2 

Example 6.8. We can cite various Lennard-Jones type potentials used in molecular simulation or 
in the study of social aggregation (see [^), besides the classical Vlj- For instance the (12 — 10) 
potential 

V(r) = ^ ^ 

V J ^12 .^10 


describes hydrogen bonds (see 24l|i. 
A (6 — 4) potential 




is also used for finding energetically favourable regions in protein binding sites (see 2^ for details). 
Lemma 6.9. Let 1 < xi < X 2 , then function r V^J (r^) is decreasing on 
and increasing on 


0 , 


02 X 2 \ '^(^2-^1) 
OlXl 


02 X 2 \ '^(^ 2 -^ 1 ) 

j 


OlXl 


Proof. The first derivative of this function is r 1 —?• — 202 X 2 ^ 2*2 1 _|_ 2aiXir 2*1 1 and 

(O'( 0>0 4 


r > 


f 02X2 \ 2 (^ 2 -^i) 
\aixij 


□ 
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Obviously, the form of potential implies that minimizer among lattices exists. Indeed, if we 
fix the area, one of the distance in the lattice cannot be too small otherwise lattice energy goes to 
infinity (see jl. Proposition 2.3] for details). 

As in our previous work j^, the following upper bound for area such that triangular lattice is the 
unique minimizer for our energy is not optimal but the best for our method. Moreover, its upper 
bound is better than we apply Cauchy’s rule (Proposition 16.31) but the method is specific for this 
kind of potential. 


) \ 

oir(x2)/ 

the unique minimizer of Eyhj, up to rotation, among lattices of area A fixed. 


, then Aa is 


Proof. We have, by proof of Theorem 16.31 for any y > 1, 


SM = iSn {»■" + !/"■') - iSA (!/■'■ + 

where gA{y) = . We compute 

9Aiy) = - «i(^2 + XI - - ai(x 2 - = y^^-^^-^Aiy) 

where UA{y) = —_ll5^ya;2+a:i-i _ + xi — l)y‘^^^~^ — ai(rE2 — xi). Moreover 


u'Air) = (x2 + xi - ^ 


(2X2 1)Q2 X 2 —XI 

^X2-Xi 9 


ai(2xi — 1) 


We have 
If A < vr 


ai(2xi — ^ /air(x2) \ /2xi —lW2-2:i 


^A^y) ^ y 02(2x2-1) ) vr V«2 r(xi)y V2X2-1, 

———^ then y <1 and u'Ay) > 0 on [l;+oo), i.e. ua is an increasing function 

air(x2)y 
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on [l;+oo). Furthermore we have 


'Wyl(l) 


( 2 X 2 - 1 ) 


0:2 

y^X2—xi 


— Oi 


(2X2 - 1) 


02 ^ 
^X 2 -xir(x 2 ) 


ai7r^2 1 

r(xi) . 


> 0 


and g'J^ is positive on [l,+oo). Thus gA is increasing on [l,+oo) and, always by assumption, 


sMl) = 2 - a,) > 0 


hence gA{y) > 0 on [l,+oo) and by Theorem 11.11 is the unique minimizer of EyLj, up to 
rotation, among Bravais lattices of fixed area A. □ 


Remark 6.11. This bound is optimal for our method because, gAi^) 
and A I— )• 51 ^( 1 ) is a decreasing function. 


0 for A = TT 


/ a 2 r(xi) \ 
Voir(x 2 )y 


1 

X2-XI 


Example 6.12. For V{r) = -^ which corresponds to Lennard-Jones energy in our case in 

[^, we hnd 

/ a 2 r(xi) \ " 2-^1 _ / F(3) _ TT 

Uir(x2) ) “ V2r(6) J ~ 1201/3 • 


Now we prove that for small parameters, the global minimizer among all Bravais lattices - without 
area constraint - of our energy is unique and triangular. We follow some ideas from our previous 
paper which cannot be apply for classical Lennard-Jones potential Vlj(x) = — 2r“®. 


Lemma 6.13. (Upper hound for global minimizer’s area) Let La^x o global minimizer of 
Eyhj among all Bravais lattices, then 

^ a,x 


|Ta,a;| ^ 


02X2 \ ^2-^1 
aixij 


Proof. Same argument of STEP 3 in the proof of Theorem II. 21 B.1. 


□ 


Thus we can prove Theorem I1.2L B.2 : 

Proof. Let La^x be a global minimizer of EyLj. We have 


h{x2) < h{xi) TT 

then by Lemma 16.131 we get 


/ a2r(xi) \ ^ 2-^1 ^ / 02X2 \ ' 
Vair(x2)y “ \aixij 

1 

/ 02 r ( xi )\ ^2-^1 


Tax — 1 -p/ \ I 

Voir(x2)y 

and by Proposition 16.lOl the minimizer among lattices of area |Ta,x| fixed is unique and triangular, 
hence the global minimizer of the energy is unique and triangular. Furthermore, let 

/(r) := EyLj[rAi] = a 2 CAi( 2 x 2 )r“ 2*2 _ oiCa^(2 xi)r“2"'i 

then we have f'{r) = — 2 a 2 X 2 CAi( 2 x 2 )r“^* 2 -i _|_ 2 aiXiCAi ( 2 xi)r“^^i“i and 


f'{r) > 0 


r > 


/ a 2 X 2 CAi ( 2 x 2 ) y (" 2 -^ 1 ) 
VaiXiCAi(2xi)y 
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a2X2CAi{‘2X2) 

aiXiCAi(2xi) 


1 

X2-XI 


Hence the minimizer of EyLj among triangular lattices is 


with \La^j; 



□ 


Remark 6.14. For an easy numerical computation of global minimizer’s area, we can use formula 
(j2.2l) to obtain 

1^ I 1 / «2a;2C(3:2)(C(^2,l/3) - C(a^2,2/3)) y2-^i 

2^3 V«i®iC(3:i)(C(a;i,l/3) - C(a;i,2/3))y 

Remark 6.15. We can apply the previous Theorem to x = (2, 3), and r (r^) is a (6 — 4) 

potential. Moreover we can choose ai and 02 such that the well is as deep as we want. 



Now we explain a method to choose xi,X 2 in order to have a triangular global minimizer and we 
give several numerical values. 

Lemma 6.16. (Variations of h) Function h is a decreasing function on [l,y;“^(log7r) — 1) and 

T'ix) 

increasing on ['i/;“^(log7r) —1, + 00 ) where i({x) = ( is the digamma function defined on (0, + 00 ). 

r(x) 

Proof. We have h'{t) = r(t) + tT{t) — t log7rr(t) and 


h'{t) > 0 


V’(i) + f ^ 


We use the famous identity ^{t) + j = fi(l-\-t) for any t > 0 and we obtain, because is increasing 
on ( 0 , + 00 ), 

h!{t) > 0 t > y^“^(log7r) — 1. 


□ 


Remark 6.17. We compute '!/;“^(logvr) — 1 ss 2.6284732 and we define M 1 such that h{M) = 
h{l). We have M ~ 4.6022909. Thus, if we want apply the previous theorem, it is clear that 
xi < '(/’“^(logvr) — 1 and X 2 < M. Moreover, if we choose xi € (l,y)“^(log 7 r) — 1), we can choose 
X 2 £ {xi,Mxfi) where 7 ^ xi is such that h{Mxfi} = h{xi). 
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1.5 


2.5 


3.5 


4.5 


tp ^(logTr) — 1 


M 


Fig. 8 : Graph of h 

Unfortunately we can only choose X 2 and xi such that 1 < xi < X 2 < 4.6022909 and Lennard Jones 
case (x 2 = 6 and xi = 3) is not covered by our Theorem II. 2I B.2. 


We compute, for different values of (xi,X 2 ) satisfying h{x 2 ) < h{xi) and for a = (1,1) 


the value of the minimizer of y i-e. := 


X2 \ 2(a:2-3:i) 

Xl 


the value of the length of triangular global minimizer, i.e. Ta,x) he. := , 

the densit}!!!] of i.e. da,x ■= \La 


1-1 


Xl 

X2 

1.1 

1.5 

2 

1.5 

(1.47,0.64,2.78) 



2 

(1.39,0.80,1.82) 

(1.33,0.95,1.27) 


2.5 

(1.34,0.90,1.45) 

(1.29,1.02,1.10) 

(1.25,1.10,0.96) 

3 

(1.30,0.95,1.27) 

(1.26,1.06,1.03) 

(1.22,1.11,0.93) 

3.5 

(1.27,0.99,1.19) 

(1.24,1.08,1.00) 


4 

(1.25,1.01,1.14) 




Table 1. Some values of {y^'^,ra^x,da 


Obviously, we have non-optimality of if A is sufficiently large, given by Proposition 16.61 : 

Proposition 6.18. (Lennard-Jones at low density) Triangular lattice Ayi is a minimizer of 
Eyhj among lattices of area A fixed if and only if 

^ f Q2(Cl(2X2) - CAi(2X2)) \ -2-1 

~ |L|=i,L^Ai Vai(CL(2xi) - CAi(2xi))y 

i.e. if A is suffieiently large, Aa is not a minimizer of EyLj among lattices of fixed area A. 

^^Here we exceptionally give values of densities - and not areas - more used in molecular simulations. 
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Proof. We apply directly Proposition 16.61 with = 1 and 


inf max 

i£l, 

\L\=1 


f tl{/+}ai(CL(2xi) - CAi(2xj)) \ i 

V I«i|(Cl(2xi) - Cai(2xi)) ) 


f Q2(Cl(2x 2) - CAi(2a:2)) \ " 2—1 
|L|=i%Ai Vai(CL(2xi) -CAi(2a;i))y 


□ 

Remark 6.19. More precisely we can found an explicit computable bound (but not optimal) if 
we take L = 7? and use m and (EJD. We give here densities do such that for any 0 < d < do, 
EyLjld-^/^Z"^] < EyLj[d-^/^Al], i.C. square lattice have less energy than triangular lattice, with 
ai = 02 = 1. 


Xl 

X2 

1.1 

1.5 

2 

3 

4 

5 

6 

7 

8 

9 

1.5 

0.05 










2 

0.14 

0.31 









2.5 

0.21 

0.37 

0.43 








3 

0.27 

0.41 

0.47 








3.5 

0.31 

0.45 

0.50 

0.58 







4 

0.35 

0.48 

0.53 

0.61 







5 

0.42 

0.53 

0.58 

0.65 

0.71 






6 

0.47 

0.58 

0.63 

0.69 

0.74 

0.78 





7 

0.52 

0.62 

0.66 

0.72 

0.77 

0.80 

0.83 




8 

0.56 

0.65 

0.69 

0.75 

0.79 

0.82 

0.84 

0.86 



9 

0.60 

0.68 

0.72 

0.77 

0.81 

0.84 

0.86 

0.88 

0.89 


10 

0.62 

0.70 

0.74 

0.79 

0.83 

0.85 

0.87 

0.89 

0.90 

0.91 


Table 2. Non-optimal critical densities for non-optimality of triangular lattice. 


7 Potentials with exponential decay 

7.1 Definition and prove of Theorem II.IL A for fa,x,b,t 

Definition 7.1. Let a = {ai,...,an) G (M*)” with an > 0, x = {xi,...,Xn) be sueh that 3/2 < xi < 
X 2 < ... < Xn, b = {bi, ..., bm) G (M*)”* and t = {ti, .■.,tm) G (M*)”*, we define 


n 

fa,x,b,t{r) ■■= 

i=l 


3=^ 


m 

We set I- := {i‘,ai < 0} and B := \bj\tj. 

i=i 

Remark 7.1. As explained in [^, Fumi and Tosi proposed in a potential for interaction 
between ions Na~^ and Cl~ defined by 


V{r) 


— + 


J.6 


03 


Obviously, potential r i-)- — is not admissible but the form of V is close to faxbt- 

r ’ ’ ’ 

Let us prove Theorem II.2I A for this potential. 
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Proposition 7.2. If we have 


A < min < vr min ■ 


anT{xi) 


driT^ 


Xn + l 


’ V(K^-} + i)5r(x„) 


3371 + 1/2 


(7.1) 


is/- |^V(2tt{l-} + 2)|ai|r(x„), 

then Aa is the unique minimizer of up to rotation, among Bravais lattices of fixed area A. 

Proof. As we have, by classical formula, for a > 0, 

a 
la 


£-l[e-^^](y) = ^y-3/2e 
zy vr 


OiTT 


Xi-l 


r(xi) 


we obtain, for any y > 0, 


taking a = tj for any 1 < j < m and setting Oj = 

n ^ t ^ ^ TD 

and it follows that 

9A{y) ■■= y~^hfa...b.t +hfa,.,b,t (^) 

. 5A3/2 5yl3/2 

^ 7 , J- \ nr. 


- 3/2 


E+Tte-'‘+: 


2=1 


27r2 27r2y3/2 


= y 


I X„+a:i-l', _ ^ Xn + 1/2 _ ^ 3:n-3/2 

1 27r2 ^ 27r2 ^ 

. 2=1 


We set 


Pa,; 


" R/13/2 R/13/2 

Mt!') E 

2 = 1 


and we notice that, for any 1 < z < n, 

^2 7^ 3/2 

Xi 1 ^ Xtt, “h 1/2 

Xn Xi 1 ^ Xtj, 3/2 

Qin .2Xr7, — 1 

la;n —1 U 

negative terms is 2(|{/_} + 2. Thus, by Cauchy’s rule (12.51) . an upper bound on the values of the 
positive zero of Pa,x,b,t is 


because Xi > 3/2. Hence the higher order term is ^ with an > 0, and the number of 


Mn max < max 

Pa,x,b,t ' 


/(2K/-} + 2)|a/A- 


^i \ xn+x.; — ! 


Q^r, 


, max 
2G7_ I 


/(2tt{/-} + 2)|a/A- 


n \ xn — XA 


Q-r 


I 2g 7_ I y 

( H(2tt{I_} + 2 )A’""+i/ 2\ /5(2tt{I_} + 2)A"’"+^/2\ 




Now we have 


A < min < tt min 


27r2Q;^ 


®nr(xj) 




iei-\ V( 2 tj{/-} + 2)|ai|r(xn)/ [ ’ V (2tl{7-} + 2)Hr(xn) 


27r2Q;,; 


2a„Tr^'^~^^ \ "“"+1/2 
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Vi & I-,A < TT 




(2tj{/_} + 2)|ai|r(x,; 


and A < 


2a„7r^"+i 


m{I.} + 2)BT{xn) 


a^n + 1/2 




^ 5(2H/-} + 2)A^"+V2 

and -- < 1 


Q!.r, 


27r^ 


Q!.r, 


^Pa,x,b,t — 


therefore if ?/ > 1 > then Pa,x,b,t{y) > 0 hence gA{y) > 0 and by Theorem 11.11 is the 

unique minimizer of Ef^ ^ ^ ^ among Bravais lattices of fixed area A. □ 


Corollary 7.3. If I- = $ and 


A < 


/ 0,1 \ a:„ + l/2 


V Br(xn 

then Aa is the unique minimizer of among Bravais lattices affixed area A. 

Remark 7.4. Obviously, for any Aq, there exists B sufficiently small such that for any A G (0, Aq], 
is the unique minimizer of our energy among Bravais lattices of fixed area A. We will study a 
simple particular case in next subsection in order to illustrate this fact. Furthermore we skipped 
the completely monotonic case but in the next following part we will give explicit condition for 
complete monotonicity in a simple case (see Proposition 17.6p . 


7.2 Example : opposite of Buckingham type potential 

In this part we study opposite of Buckingham type potential. Indeed, we cannot study Buckingham 
potential 

Tr / \ —ctr ^2 ^3 

VB{r) = aie - ^ ^ 

because lim^^o yBif) = —oo and limr^+oo Vsir) = 0 and it is sufficient to do ||u|| — >• 0 in order to 
have Evg[L] —)■ —oo. Hence we choose to treat simple general approximation of its opposite, well- 
adapted for our problem of minimization among Bravais lattices. Moreover we simplify notations 
in order to have only two parameters : 

Definition 7.2. For a = ( 01 , 02 ) G (0,-|-oo)^ and for x = (xi,X 2 ) G (0,-|-oo) x (3/2,-boo), we 
define, for r > 0, 

/a,x(r) = 02r-"^-oie-"i^^. 


Lemma 7.5. (Variations of potential r t-A /a,x(?’^)/ fke have the following two cases : 
1. ^/( 2 x 2 + 1) 


In I +i±i I - 1 


S. if (2X2 + 1) 


In 


Tm < 


2x2+1 

xi 


Xi 

2X2 + 1 

Xi 


< In ( then r i-+ faxi'i^'^) is decreasing on (0,+oo); 

V / 


— 1 > In (—^ ) then there exists rm,rM £ (0,+oo) such that 

V / 

< Cm ond r i-+ fa,x{'f‘^) is decreasing on intervals (0,rm) and {rM,+oo) and 


increasing on {rm,rM)- 

Proof. We have /(r) := fa,x{r‘^) = a 2 r~‘^^^ — aie~^^'" and 


pf / \ 2O2X2 —XTT 

fir) = -3^ + aia;ie 


Thus we get 


fir) > 0 


.^2X2 + 1 


-afir^2a;2+l ^ 


202 x 2 

OlXl 


gir) > 0 
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where 


g{r) = —xir + ( 2 x 2 + 1) In r — In 


202X2 


As g' (r) =-, g is increasing on ( 0, j and decreasing on ( , +oo j. Moreover 

g{r) goes to —oo as r —)■ 0 or r —)■ +oo. 

Hence if g < 0, i.e. 


(2x2 + 1) In 


2X2 + 1 


— 1 < In 


202X2 


then g{r) < 0 and f'{r) < 0 on (0,+oo), i.e. / is decreasing on (0,+oo). 

Furthermore, if g (> 0 then there exists rm,rM such that < vm and / is 

\ X\ J XI 

decreasing on intervals (0, r^) and {rM,+oo) and increasing on {rm,rM)- 


□ 



Fig. 9 : Graph of r i-)> fa,x{r‘^) for a = (1,1), x = (5, 6) on the left and a = (1,2), x = (1,6) on 

the right. 


Proposition 7.6. We have the following two eases : 


• If it holds 



> In 


/ OiXiF(x2) \ 

V 2^/^:a2 J 


then for any A > 0, is the unique minimizer of up to rotation, among Bravais 

lattiees of fixed area A. 



< In 


/ OiXiF(x2) \ 

V 2^/^:a2 ) 


then Ayi is the unique minimizer of Ef^^, up to rotation, among Bravais lattices with fixed 
area A. Moreover, for any a G (0, +oo)^, X 2 > 3/2, Aq > 0 and any xi such that 


0 < xi < Cao ■= 


02^^+^ 

OiAq^’^^/^F(x2) 


A/i is the unique minimizer of Ef^^^, up to rotation, among Bravais lattices of fixed area 
A G (0, Aq]. 
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Proof. By classical formula, we get 


0-2 


™2 

-— 1 . 

2,2-1 _ aiXi ^_y2^ Ay 




Our theorem is a consequence of Proposition 13.11 because 


Vy > o>/^/a..(y) > 0 (^2 +1/2) 

Indeed, we have 


1 + In 


™2 

Xi 


4:X2 + 2 


> 


V 2^/^^a2 J 


Vy > 0,t^fUy) > 0 ^ > 


.a:?/4j/„a;2+l/2 > QlXir(x2) 
2-/7ra2 


We set 


<,(!/) = # + fe + l/2) In, 


Ay 


V 2^/^Ta2 J 


and we have g'{y) = ■ It follows that <7 is decreasing on ^ 0 , increasing on 

( 4 x 2+2 ’ + 0 °) • 9 Soes to +00 as y goes to 0 or + 00 , it is clear that 


Vy > 0,y(y) > 0 


4x2 + 2 


{X2 + 1/2) 


> 0 


1 + In 


4x2 + 2 


> In 


/ aixir(x2) \ 

V 20ra2 j 


Now, if fa,x is not completely monotonic, we apply directly Proposition 17.21 to obtain second point. 
Third point is clear because for any ( 01 , 02 ) E (0,+oo)^ and any X 2 > 3/2, 

/ 027r*^'''^ \ ^ 2 + 1/2 

VoiXir(x2) J 


is an increasing function which goes to infinity as xi —>■ 0 . □ 

Example 7.7. For instance, we can choose a = (1, 1), X 2 = 6 and ^0 = 1- Thus we get 

Cl = — « 0.0727432 
11 ! 

and for any xi < Ci, Ai is the unique minimizer of among Bravais lattices of unit fixed area. 
Moreover the form of the potential y fa,x{y‘^) is such that the decay to 0 at infinity is slow as xi 
goes to 0 . 
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Remark 7.8. Our argument used in proofs of Theorem 11.21 based on variations of potential, can’t 
be applied for our potentials fa^x- 

Acknowledgements: I am grateful to Etienne Sandier, Florian Theil, Salvatore Torquato and my 
colleague Peng Zhang for their interest and helpful discussions. 
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